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Abstract. Our report devotes a 3D numerical model of the interaction of the 
solar wind with the Solar Probe Plus (SPP) spacecraft. The SPP model includes 3 
main parts, namely, a non-conducting heat shield, a support system, and cylindrical 
section or spacecraft bus that contains the particle analysis devices and antenna. The 
simulation was performed using 3D hybrid code which describes the proton dynamics by 
the particle motions whereas the electrons are considered in a fluid approximation. One 
observes an excitation of the low frequency Alfven and whistler type wave directed by 
the magnetic held with an amplitude about of (0.06-0.6) V/m. The compression waves 
and the jumps in an electric field with an amplitude about of (0.15-0.7) V/m were also 
observed. The wave amplitudes are comparable to or greater than previously estimated 
max wave amplitudes that SPP is expected to measure. The results of our hybrid 
simulation will be useful for understanding the plasma environment near the Solar Probe 
Plus spacecraft. Future simulation will take into account the charging of the spacecraft, 
the charge separation effects, an outgassing from heat shield, a photoionization and an 
electron impact ionization effects near the spacecraft. We also need to perform a longer 
simulation in order to receive a ’’steady-state” solution with a formation of a plasma 
wake in quasi-parallel interaction. 

Numerical simulation studies; Key words: Solar wind; Alfven waves, Whistlers, 
Ionospheres; Atmospheres; Induced magnetospheres; Magnetic barrier; Satellites 
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1. Introduction 

The numerical simulation of the interaction of space plasma with spacecraft devices 
plays a key-role in their design. The study of the excitation of waves near a rapidly 
moving body has long history starting from the quasi-analytical approach (see e.g. 
Alpert (1974) and the references thereby) to the recent electrostatic modeling. While 
the electrostatic studies produced very important information concerning the charging 
of spacecraft, plasma void and wake, the problem of excitation of the electromagnetic 
waves upstream and downstream flow near spacecraft are great interest for onboard 
measurement in the SW, magnetosphere and ionosphere. This is especially true for SPP 
where new regimes of solar wind plasma are being explored. The simulation provides 
the general characteristics of the plasma environment and the electromagnetic field 
distribution near the probe and at it’s surface. The external surface of the ”SPP” 
consists of a conducting part and an insulating part. It moves in the supersonic/subsonic 
and superalfvenic/subalfvenic solar wind flow. 

In the first approximation the plasma environment near the ” Solar Probe+” may be 
similar to the plasma environment near the Moon with a weakly conducting surface. In 
the case of a non-conducting model of the Moon the solar wind particles penetrate the 
surface on the day-side of the Moon whereas on the night-side a plasma wake with the 
low-density void is formed. On the day side of the Moon’s surface the plasma particles 
are absorbed, and the perturbation region forms a thin boundary layer of thickness 
S ~ c/a;p e oc 1km at 1AU. (Neugebauer, 1960). On the night side the perturbation 
region is bounded by a surface of weak perturbations forming a cone with half-apex 
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angle sin 2 9 = 1 / M\ + 1/M| oc 8 and Ms ex 8 are the Alfven and sound Mach 
numbers, Va ~ 50km/s, w t h ~ 50km/s, [/ ~ 400km/s, and d ~ 10°) (Whang, 1969; 
Lipatov, 1976). 

No bow shock or oblique Mach cone are formed. The perturbation of the magnetic 
held inside the wake is determined by drift currents [Wang, 1968; Wang, 1969; Wang 
and Ness, 1970; Lipatov 1976, Lipatov, 2002], Kinetic instabilities may also play a role 
in wake dynamics. Hybrid simulation of the interaction between supersonic plasma 
how and a weakly conducting body shows formation of a strong Mach cone outside the 
plasma wake [Lipatov, Motschmann, et al., 2005]. 

2. Formulation of the Problem and Mathematical Model 

2.1. Solar Wind Simulation Model 

The interaction of solar wind particles with the ”SPP” is more complicated 
because the gyroradius of the protons (for Maxwellian core of the velocity distribution 
50 — 500 km) is much larger than the size of the spacecraft (2-3m). For energetic 
component the ion gyroradius may be extremely large with respect to the spacecraft. So 
the kinetic approximation for ions is an essential part of the mathematical model. The 
electron gyroradius may vary in range from 2.5 x 10 2 m at Earth orbit to much smaller 
scale near Sun. Therefore, electrons need a fluid-kinetic approximation to study the 
plasma environment along the ’’Solar Probe+” trajectory. 

We expect to study the plasma environment near the ” SPP” by various simulations; 
e.g., (a) standard hybrid simulation (ion in kinetic approximation, electron in fluid 
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approach) on the large scale; (b) fully kinetic implicit simulation with kinetic model 
for electrons and ions (see e.g. Hewett and Langdon, 1987; Lipatov and Lobachev, 
1996; Lipatov, 2002; Damiano, Sydora, Samson, 2003; Lipatov and Rankin, 2008) 
incorporated in the large scale hybrid model. The last simulation will take into account 
the charge separation near the surface of the spacecraft and finite electron gyroradius 
effects. We will take into account the realistic distribution of the spacecraft surface’s 
conductivity. Our simulation will serve as an expert system for design of the ’’Solar 
Probe Plus” spacecraft. The present model of the interaction of the solar wind with the 
SPP does not take into account several effects in plasma environment near spacecraft. 
Future simulation will take into account the charging of the spacecraft, the charge 
separation effects, an outgassing from heat shield, a photoionization and an electron 
impact ionization effects near the spacecraft. 

To study the interaction of the solar wind with the SPP spacecraft we use a 
quasi-neutral hybrid model for ions and electrons. 

In our coordinate system the X axis is parallel to the solar wind velocity - U 0 , Z is 
aligned with the equatorial plane, and Y = Z x X is perpendicular to equatorial plane 
and anti- aligned with the E = — x B electric held (electric held points upward). 

In the hybrid simulations described here, the dynamics of upstream ions is 
described in a kinetic approach, while the dynamics of the electrons is described in a 
hydrodynamical approximation. 

The single particle ion distribution function /)(f,x, v) has to fulfil the 
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Vlasov/Boltzmann equation 


— f-L —f _l F 9 

dt fi + V dx + 



( 1 ) 


where F symbolizes forces due to electric and magnetic fields acting on the ions, F co u 
is the collision term. In this paper we use the particle-mesh model for ion dynamics 
instead of the Vlasov equation, Eq. ( 1). 

The single ion particle motion is described by the equations (see, e.g. Eqs. (1) and 
(14) from [Mankofsky, Sudan and Denavit, 1987]): 


dr it i dvjj e / v,,i x B \ m e u ie 

dt 1,1 ) dt c ) Mieni 


( 2 ) 


Here we assume that the charge state is Z t — 1 and that all ions have the same mass 
Mi. U i and J denote the charge-averaged velocity of all (incoming and pickup) ions and 
the total current, Eq. (6). The subscript i denotes the ion population and the index l 
is the macro-particle index. v ie is collision frequencies between ions and electrons, that 
may include Coulomb collisions and collisions due to particle-wave interaction. Note 
that the collision rates used in Eq. (2) must depend on individual velocities of ions and 
electrons. However, we use the effective resistivity rj, rj = a = m e /(ne 2 T e ), where 
r e = . The electrical conductivities may be estimated as 


<7 _l = di T e 3/2 , (7 1| = 1.92erj_, a x = 0.9 x 10 13 /((A/10)Z;) s" 1 • eV" 3/2 , (3) 


where T e denotes the electron temperature and A is the Coulomb logarithm (see, 
e.g., pages 215-216 from [Braginskii, 1965]). For the typical solar wind parameters 
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T e = 100 eV (electron temperature) and no = 10 4 cm -3 (density) the electrical 
conductivities are cr± ~ 4.7 x 10 13 s^ 1 and cry ~ 9.2 x 10 13 s _1 . 

In our simulation we use the low (much smaller than the real value) effective 
conductivity to suppress ’’shot” noise and for modeling Solar Probe Plus’s bus; hence, 
we must keep the first collision term in the right hand side of Eq. (2). 

In the nonradiative limit Ampere’s law is given by 

4:71 

— J = V x B; (4) 

c 

and the induction equation (Faraday’s law) by 

+ VxE = 0. ( 5 ) 

c at 

The total current is given by 

J J e + Jj, J) (6) 

We further assume quasi- neutrality 

n e « rq. (7) 


The equation of motion of the electron fluid takes the form of standard generalized 
Ohm’s law (e.g. Braginskii, 1965): 


^ 1 ^ 1 „ m e J 

E = (J c x B) Vp e H W* — 

en e c en e e ne 


m d 
e dt 


U e , 


( 8 ) 


where p e = nm e {v ®) / 3 = n e ksT e , and v' e are the scalar electron pressure and the thermal 
velocity of electrons, and the electron current is estimated from Eq. (6). The last term 



on the right side of Eq. (8) is the electron inertial term. At 5 R s we have no = 10 4 cm -3 , 
v' e = 5000 km/s, and from Eq. (3) one gets cr e i = 0.4 x 10 14 s -1 for background plasma 
parameters. 

Since we suppose that electron heating due to collisions with ions is very small, 
the electron fluid is considered adiabatic. For simplicity we assume that the electron 
pressure may be represented as: 

p e oc nf /3 . (9) 

We also assume that n e ~ rq. Otherwise, we have to calculate the electron pressure from 
heat balance for electrons (see, e.g., Braginski, 1965) taking into account the heat fluxes. 
The ion kinetic approach allows us to take into account the effects of anisotropy of ion 
pressure, the finite gyroradius effects, the asymmetry of plasma flow around the Solar 
Probe Plus, and to estimate the particle flux at the surface of the spacecraft. Remember 
that the fluid models which account only for the scalar (i.e., isotropic) ion pressure may 
result in an extra-expansion of the ions along the magnetic field. Our plasma model 
may support the kinetic ion cyclotron wave and electron fluid cyclotron waves. 

2.2. SPP Spacecraft Model 

Figure 1 shows a scheme of the SPP spacecraft, a system of coordinates and a 
direction of the plasma flow and magnetic field. The SPP model includes 3 main parts, 
namely, a non-conducting heat shield (1), a support system (2), and cylindrical section 
or spacecraft bus (3) that contains the particle analysis devices and antenna. The heat 
shield has the following geometrical parameters: the diameter, Dg^id = 2.7 m, and the 
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thickness = 0.335 m. The gap between the heat shield and cylindrical section is equal 
1.188 m. The cylindrical section or bus has a diameter equals 1.026 m and a length 
1.188 m. We also take into account the effective resistivity of the Solar Probe+’s bus, 
the heat shield and the gap between the heat shield and the spacecraft bus: 

Pbus = (1 — 15) x 10 -3 ohm-m-, p s hieid = (3 — 15) ohm-m] p gap = (3 — 15) x 10 ~ 3 ohm-m 

( 10 ) 

Our code solves equations (1) - (5), (6) - (9). Here we note that the gap has several 
trusses to provide mechanical interface between bus and heat shield. If electrical 
conductivity of trusses high enough that differential charging between spacecraft bus 
and heat shield is low then the heat conduction from heat shield to spacecraft bus 
maybe too high. We are not aware of the exact number at this time with regard to truss 
electrical conductivity parameters. 

Initially the computational domain contains only supersonic (subsonic) solar wind 
flow with a homogeneous spatial distribution and a Maxwellian velocity distribution. 
The magnetic field and electric fields are B = Bo and E = Eo = — Uo x Bo- Inside Solar 
Probe spacecraft electromagnetic fields are E =0 and B = B 0 , and the bulk velocities of 
ions and electrons are also equal to zero. In the cases examined here we choose spiral 
angle of field 9b u = 11° — 45°. 

Far upstream (x = —DX/2), the ion flux is assumed to have a Maxwellian 
distribution, 

/ = exp 


( v ~ U ) 2 
2 < 


( 11 ) 
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where Vth and U are the thermal and the bulk velocities of the solar wind plasma flow. 
We have not included the importance of the field aligned strahl electrons which are 
non-Maxwellian. 

Far downstream, we use a buffer zone to provide the return of the particle with 
negative value of the u velocity component. We also adopted Sommcrfcld’s radiation 
condition for the magnetic held. On the side boundaries (y = ±DY/2 and z = ADZ/ 2), 
unperturbed boundary conditions were imposed for incoming how particles and the 
electromagnetic held. At Solar Probe Plus’s surface the particles are rehected or 
absorbed. There is no boundary condition for electromagnetic held, and we also use 
our equations for the electromagnetic held, Eqs. (4-5) and 8 inside Solar Probe Plus 
spacecraft but with internal conductivity and the bulk velocity that is calculated from 
the particles. In this way the jump in the electric held is due to the variation of the 
value of the conductivity and bulk velocity across the Solar Probe Plus bus’s surface. 
Note that the position of the center of the bottom of the heat shield of Solar Probe Plus 
is x = 0, y = 0, z = 0. 

The three-dimensional computational domain has dimensions DX = 10L, 

DY = 8 L, and DZ = 8 L, where L = D s hi e id/2 = 1.35 m is a radius of heat shield. 

We used mesh of 401 x 301 x 301 grid points, and 9 x 10 8 particles for protons for a 
homogeneous mesh computation. The time step for the particle update, A t p satishes the 
condition v max A t p < min (Ax, Ay, Az)/16, whereas the time step for the electromagnetic 
field time integration Atf satisfies the condition v max Atf < min(Ax, Ay, Az)/6400. 
Note that the grid spacing are the following: Ax = Ay = A z = 0.036 m. 
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The relationship between dimensional (U, E, B, p e , n, T) and dimensionless ([/', 
E' , B' , p' e , n' , T') parameters may be expressed via the dimensional upstream values as 
follows: 

U = UU, E = E'BqUo/c, B = B'Bq, Pe = Pe'Pe 0, 

n = rino, T = T'M-JJq, (12) 

whereas the dimensional time and distance may be expressed via the bulk velocity Uq 
and characteristic scale L = D shie i d / 2: 

t = t'L/Uo, x = x'L. (13) 

The global physics in SPP’s environment is controlled by a set of dimensionless 
independent parameters such as Alfven Mach number M\, the ion and electron plasma 
betas P e , the electron/proton mass ratio m/M p , diffusion lengths, and the ion 
gyroradius e = p c i/L. Here p c \ = Uo/(eB /M- X c) = M^cju) v \ and the ion plasma frequency 
ca P i = Airrioe 2 / AE The actual value of the proton gyroradius is about (0.31 — 2.5) km 
using the above formulas using the above formulas. The grid spacing has the value 
A x = L/50 = 0.036 m. In order to study the ion kinetic effect, the ion gyroradius has to 
be resolved on the grid. 

2.3. Numerical Method 

We employed a standard particle-in-cell (PIC) method with a homogeneous grid. 
The time integration of the particle equations of motion uses a leapfrog scheme. The 
time integration of the electromagnetic field equations uses an implicit finite difference 
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scheme (see, e.g., Lipatov (2002)). 

The finite difference approximation for the electric field and the electron velocity 
equations may be produced by the following way: 


dVx(Vx E n+e ) + AE n+0 + (V x E n+0 ) x I 


n+9 


; \n+6\ 


+ 0 ((V X (V x E n+0 )) x B n ) = Q, 


where 


, me 2 . . , , 9 Ate 

d = — + 9n e l d At , A = n e , g = 


M Ml 


M\ ’ 


_ OeAt { Ml lA n+ i /2 „ 


x B r 


N s 


Q = - 


o ^ t — t t \ n. \ t ~%n. tfte 


£ - W v x B ”j x B " ' + "' (;v x B ' 


and 


U. = Ui - «J/(M|n e ). 


(14) 

(15) 

(16) 

(17) 

(18) 


In the case with adiabatic electrons one can split the total electric field into the 
sum of inductive (E\ ) and electrostatic (E 2 ) fields: 


E — Ei + E 2 , 


(19) 


where E 2 satisfies the condition 


V x E 2 = 0. 


( 20 ) 


Then we can solve (14) for component E l neglecting Vp e - The electrostatic electric field 
E 2 can be calculated from p e (Lipatov, 2002) because 


Eo = 


en P 


" V p e - 


( 21 ) 
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Let us consider two meshes. The first mesh contains the nodes i,j,k, which are 
located in the center of a cell. It is used for the computation of the density, current, 
bulk velocity, electron pressure and inductive electric held. The second mesh contains 
the nodes i ± 1/2, j ± 1/2, k ± 1/2, which are located at the corners of a cell. This 
mesh is used for the computation of the magnetic held components, electrostatic held 
and final electric held. Let us assume that Jf' +1,/2 , u ve and rq are known at time level 
n + 1/2. Then, expressing V x E™ +e and V x (V x E ™ +9 ) via central finite differences 
at each cell center, one can obtain the following 3x3 matrix equation: 

L ■ gfi* = F (g£fj ±lik±t ) . (22) 

Equation (18) may be solved by iteration. In each iteration the electric held on the right 
side is given. The iteration continues until some convergence criterion is satished. At 
the same time the electrostatic held E -2 is calculated at from Eq. 21. 

The second term on right side of Eq. 17 keeps a large parameter, so that a strong 
shot noise in a density computation may cause a strong oscillation in the electric 
held. Note that the hrst term (electron inertia) on the left side of the Eq. 14 has a 
large coefficient so that this term may play very important role in damping of the low 
frequency waves. In this report we assume that m/M p = 0 and we neglect the electron 
inertia term (hrst term on the left side) in Eq. 14. 

We used different time steps for particle and held pushing (subcycling). The code 
was optimized for parallel computation using OpenMP parallel environment. However, 
the staff from the Computational and Information Sciences and Technology Office 
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(NASA GSFC) is working on the conversion the our code to the MPI environment. 
That may allow us to use a distributed memory parallel computer system for our 
computation. 

Since the gyroradius must be resolved, a grid point spacing of less than 1 gyroradius 
is required in order to avoid numerical dispersion and dissipation. On the other hand, 
good statistics are required, therefore a sufficiently large number of particles per cell 
have to be used (i.e., to obtain low “shot” noise, which manifests itself as fluctuations 
in the numerical plasma parameters due to a small number of particles per cell). 

A multiscale simulation will be conducted using adaptive mesh and particle 
refinement techniques with composite grids. This code exploits a splitting-of-particles 
procedure to keep a low level of “shot” noise on the finer mesh. At the end of the global 
time step we use the synchronization for the values of the electromagnetic held on coarse 
and finer meshes. A.S. Lipatov has developed a version of the Complex Particle Kinetic 
method (see, Hewett, 2003) [in collaboration with D.W. Hewett (LLNL) and M.R. 
Combi (Univ. of Michigan)] for the case of multiple plasma beams (inter-penetrating 
how) which allows us a reduction in computational resources by factors of 100 - 1000 in 
comparison with standard PIC simulation (Lipatov, 2008a). These codes were optimized 
for parallel computation using MPI and OMP. 

3. Results of the Simulation 

To study the interaction of the solar wind with the SP+ we use the following 


sets of the solar wind plasma parameters: = 0.1; /3 e = 0.000125 — 0.00125 — 0.1; 
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magnetic field, B 0 = 1500 — 6250 nT; bulk velocity U 0 = 400 — 200km/s, density 
n sw = 5 x 10 3 — 10 4 cm~ 3 , Alfven Mach number Ma = 1.2 — 10 and = 11° — 45°. 
Table 2 shows the all parameters that correspond the case: (a), (b) and (c). These 
parameters correspond to the values of the motional electric held Eq = UqBq « 1.24 V/m 
for r = 5 R s and E 0 0.3 V/m for r = 9.5 R s . The higher values of the magnetic held 
and density correspond to the distance from the Sun of about 5 R Sl whereas the lower 
values correspond to the distance of 9.5 R s . We present here the results of simulation 
for the first case. The time step for electromagnetic held update is 48 times smaller in a 
quasi-parallel case than it the time step in an oblique case. In the present simulation, we 
use a reduced value for the electron plasma beta, /3 e = 0.001 to reduced the oscillation 
in the polarization electric held because of the large value of the dimensional proton 
gyroradius, e = p p / L = 1840 — 230, see the equations (13-17). 

3.1. Quasi-parallel Interaction of the Solar Wind with Solar Probe+ 

Let us consider hrst the global picture of the interaction of the solar wind with 
Solar Probe-|- in the quasi-parallel case, 6b u = 11°. We have performed the simulation 
of this case only for a short time because of the small time step that needs to be done 
for numerical stability. One simulation takes around 1 month for computing with 16 
processors on the ” shared memory” Palm/Explore machine. Figure 2 shows 2D cut 
for proton density in the y — x and z — x planes. We see the picture of the solar 
wind how around the spacecraft. Behind the heat shield a density prohle forms a cone 
dne to expansion of the external plasma into the gap between the heat shield and the 
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spacecraft bus. Behind the spacecraft bus there is no lunar-type plasma wake due to 
low value of the Mach number, Ma = 1.2 — 1.5, and high values of the plasma betas, 

/ 3 P = 0.1. f3 e = 0.000125 (case a, Table 2). Possibly, the plasma wake may form in longer 
simulations. 

In this case the current closure near a spacecraft is very complicated and is directed 
by the external magnetic held. Figures 3, 5 and 7 show the 2D cuts for the electric held 
at t — 0.014T ce . Note that T ce is the electron gyroperiod. In upstream one sees the 
linear perturbation in the electric held in form of whistler/Alfven waves in upstream 
and downstream regions at the beginning of simulation t = 0.014T ce . The value of 
perturbation in the electric held is about of 5E ps 0.02 E 0 , E 0 ps 1.24V/m). 

Figures 4, 6 and 8 show 1-D cuts for the electric held prohle along the x,y and z 
axis through the point x = 1.5L, y — 0 and z — 0. 

Far from the spacecraft, the value of perturbation in the electric held component, 
E x , is about of 5E X « (0.001 — 0.005)i7o in the y- direction, whereas 5E X « 0.01-Fo 
in the z- direction. The value of perturbation in the electric held components, E y 
and E Zl is about 5E y « 0.5i7o (whistler type waves) in the x- direction (upstream 
and downstream), and 5E y « O.Oli^o in the y- direction, and 8E y « O.Oli^o in the z- 
direction. The value of perturbation in the electric held component, E z , is about of 
5E Z ~ 0.05i7o in the x- direction, SE Z ps 0.05i7o in the y- direction, and 8E Z ps 0.1i7o in 
the z- direction. 

Near spacecraft, the E x prohle has a jump, 5E X = — 1000i7o behind the heat shield. 
Near spacecraft one sees an excitation of strong Alfven waves with an amplitude , 
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8E X ~ 0.005-Eo along y axis, and 5E X ~ (0.02 — 0.15)i?o along the z axis. We also see a 
jump in E y component of the electric field, 5E y & 0.1-F 0 i n V direction and SE y ~ 0.05_E 0 
in z direction. The components E y and E z have a jump, SE « (0.02 — 0.5) £70 in !J and z 
directions. 

Figures 9, 10 and 11 demonstrate 2D cuts for the magnetic field. At the front 
of the heat shield a formation of the magnetic held barrier is observed. At the side 
parts of the computational domain compression waves were observed in simulations. 
The value of perturbation in the magnetic held is about of SB « 0.02 B 0 (« 120 nT for 
B 0 ~ 6200 nT.). 

Let us consider the result of simulation at the later time, t = 0.29T ce 
(t = 0.04T( mns j t ), where Transit is the time for particle transit from the left boundary to 
the right boundary of the computational domain. Figure 12 shows a 2D cuts for proton 
density in the y — x and z — x planes. The perturbation in the electromagnetic held 
reaches the saturation level, since the density distribution showed no change during the 
simulation. 

Figure 13, 15 and 17 show the 2D cuts for the electric held at later time 
t = 0.29 T ce . Since the density distribution had not changed during the simulation and 
the perturbation in the electromagnetic held reaches the saturation level (left hanging). 
We do not see any strong transverse perturbation in the electric helds in the upstream 
region, but we do see large perturbations in the downstream region. The value of these 
perturbations is about of SE ps 1.0 A 0 « 1.24V/m. 

Figures 14, 16 and 18 shows 1-D cuts for the electric held prohle along the x, y and 
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z axis through the point x = 1.5L, y = 0 and z — 0. 

Far from the spacecraft in the downstream region, the value of perturbation in the 
electric held component, E x , is about of SE X ~ 0.0055 0 in the y- direction, whereas 
SE X « 0.0 l^o hr the z- direction. The value of perturbation in the electric held 
components, E y and E ZJ is about of SE « 2.0E 0 (whistler type waves) in the x- direction 
(downstream), Fig. 14, and 5E y « 0.055o in the y- direction, and SE y « 0.035o hi the 
z- direction, Fig. 16. The value of perturbation in the electric held component, E z , 
is about of 5E Z « 0.055 0 in the x- direction, SE Z ~ 0.0255 0 in the y- direction, and 
5E Z « 0.0255 0 in the z- direction, Fig. 18. 

Near spacecraft, the E x prohle has a jump, SE X = 2505 0 behind the cylindrical 
section (conducting bus). We also see a jump in E y component of the electric held, 
SE y ~ (0.5 — 0.75)Eo in y direction and SE y ~ 0.2Eq in z direction. The component E z 
has a jump, SE rs (0.15 — 0.2)E 0 in y and z directions. 

Figures 19, 20 and 21 demonstrate 2D cuts for the magnetic held. At time, 
t = 0.29T ce (t = 0.04 T transit ), we do not observe any magnetic held barrier near the 
heat shield. At the side parts of the computational domain compression waves were 
observed in the simulations. The value of perturbation in the magnetic held is about of 
SB « 0. 0250- 

Let us consider the results of the simulation with higher value of f3 e = 0.025 (case 
b, Table 2). Figure 22 shows a 2D cuts for proton density in the y — x and z — x 
planes. We see the picture of the solar wind how around the spacecraft. Behind the 
heat shield a density prohle forms a cone due to expansion of the external plasma into 
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the gap between the heat shield and the spacecraft bus. Behind the spacecraft there is 
no lunar- type plasma wake due to low value of the Mach number, Ma = 1.2 — 1.5. 

The general picture of the interaction of the solar wind with spacecraft is very 
close to the picture that was described above. Figure 23, 25 and 27 show the 2D cuts 
for the electric held at t — 0.014T ce . Figures 24, 26 and 28 shows 1-D cuts for the 
electric held prohle along the x,y and z axis through the point x = 1.5 L,y = 0 and 
z = 0. In upstream and downstream we observe the linear perturbation in the electric 
hied - whistler/ Alfv’en waves. The value of perturbation in the electric held (E y ,E z ) 
in upstream and downstream is about of SE « 0.477o (whistler type waves) in the re- 
direction, and the wavelength much larger than in case with f3 e = 0.000125. 

Far from the spacecraft, the value of perturbation in the electric held component, 
E x , is about of 5E X ps 0.002f7o in the y- direction, whereas 8E X ps 0. 005770 in the re- 
direction. In other directions we observe perturbation - 5E y ps O.Ol-Eo in the y- direction, 
and SE y ps 0.0177o hi the z- direction. The value of perturbation in the electric held 
component, E z , is about of SE Z « 0.05-E 0 in the x- direction, 5E Z « 0.01f7 0 hi the y- 
direction, and 8E Z ps 0.0 IAq in the z- direction. 

Near spacecraft, the E x prohle has a jump, 5E X = —20Eq before the heat shield 
and 5E X = 30i7o behind the cylindrical section. Near spacecraft we see an excitation of 
strong oblique Alfven (magnetosonic) waves with an amplitude , SE X ps (0.005 — 0.01) -Co 
along y axis, and 8E X ps (0.02 — 0.15)A 0 along the z axis. We also see a jump in E y 
component of the electric held, 5E y & 20 E 0 in y direction and SE y fa 0.05E 0 in z 
direction. The components E z has a jump, 5E ps 0.04i7o in the y and 8E ps 20Eq in the 
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z directions. 

Figures 29, 30 and 31 demonstrate 2D cuts for a magnetic field. At the front of the 
screen a formation of a magnetic field barrier was observed. At the side parts of the 
computational domain compression waves were observed in simulations. The value of 
perturbation in the magnetic field is about of 5B « 0.02B o . 

Let us consider now the result of simulation with (3 e = 0.025 at the later time, 
t = 0.03T ce (t = 0.01 T transit ), where T transit is the time for particle transit from the 
left boundary to the right boundary of the computational domain. Figure 32 shows 
a 2D cuts for proton density in the y — x and z — x planes. The perturbation in the 
electromagnetic field reaches the saturation level, since the density distribution has not 
changed during the simulation. 

Figure 33, 35 and 37 show the 2D cuts for the electric field at t = 0.03T ce . Since 
the density distribution had not changed during the simulation the perturbation in the 
electromagnetic field reached the saturation levels. Figures 34, 36 and 38 shows 1-D cuts 
for the electric held profile along the x,y and z axis through the point x = 1.5 L,y = 0 
and z — 0. 

Far from the spacecraft, the value of perturbation in the electric held component, 
E x , is about of SE X « (0.01 — 0.015)i7o in the y- direction, whereas 5E X & O.Ol-Fo in 
the z- direction. The value of perturbation in the electric held components, E y and 
E z . is about of 5E y « 0.02i7o (whistler/Alfv’en waves) in the x- direction (upstream 
and downstream), and 5E y « 0.01-Fo in the y- direction, and SE y « 0.01-Eo in the z- 
direction. The value of perturbation in the electric held component, E z , is about of 
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5E Z r* 0.05i?o in the x- direction, 5E Z r; 0.05-Eo in the y- direction, and 5E Z r* 0.1 Eo in 
the 2 - direction. 

Near spacecraft, the prohle has a jump, 8E X = — lOOOEo behind the heat shield. 
Near spacecraft we see an excitation of a strong Alfven waves with an amplitude , 

5E X R! 0. 005E 0 along y axis, and 8E X fa (0.02 — 0.15)E 0 along the z axis. We also see a 
jump in E y component of the electric held, SE y ~ 20Eq in y direction and SE y ~ 0.05Eo 
in z direction. The component E z has a jump, SE rs (0.02 — 20.)E , 0 in y and z directions. 

Figures 39, 40 and 41 demonstrate 2D cuts for the magnetic held. At the front 
of the heat shield a formation of the magnetic held barrier or build up was observed. 
At the side parts of the computational domain a compression waves were observed in 
simulations. The value of perturbation in the magnetic held is about of 5B ss 0.02I?o- 

3.2. Oblique Interaction of the Solar Wind with Solar Probe Plus 

Let us consider the interaction of the solar wind with Solar Probe Plus in the 
oblique case, 6f m = 45°, and Ma = 2 at time, t = 1.875 Transit, where T transit is the time 
for particle transit from the left boundary to the right boundary of the computational 
domain (case c, Table 2). The simulation demonstrates the formation of the plasma 
wake. A strong perturbation in the plasma density along the wake was observed, 

Fig. 42. The density prohle is a little bit disturbed near the side boundaries. One may 
see that SP+ bus is located inside the low plasma density cavern. This means particle 
instrument in wake may see nothing. Figures 43-45 and 46-45 show the distribution of 
the electric and magnetic helds, respectively. The asymmetry of the distributions in E 
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and B appears to be caused by the finite gyroradius effects of incoming protons. A weak 
perturbation of the magnetic field was observed near the SP+ spacecraft: compression 
of the upstream magnetic field and decompression in the plasma wake, Fig. 48. 

One can see the asymmetry in the distribution relative to the x-axis due to the 
angle between the bulk velocity and the upstream magnetic field and to the y-axis 
due to effects of the finite ion gyroradius, Fig. 47. The value of perturbation in the 
magnetic field is about of 5B « 0.02B 0) Fig. 48. Table 3 accumulates the values of the 
perturbation in the electric and magnetic field and the jumps in the normal component 
of the electric field at the surface of the spacecraft bus. Note that strong jumps at the 
bus surface may be due to the specific boundary condition and the future computations 
may reduce these values. 

4. Conclusions 

3D hybrid simulations of the interaction of the solar wind plasma with Solar Probe 
Plus, have demonstrated several new features: 

• In the quasi-parallel case the current closure near a spacecraft is very complicated 
and is directed by the external magnetic field. At the front of the heat shield a 
formation of the magnetic field barrier was observed whereas a formation of the 
strong whistler/Alfven waves were observed in both upstream and downstream 
regions. At the side parts of the computational domain compression waves were 


observed in simulations. The values of the electric field oscillation near the 
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spacecraft bus may be the same order as the maximum of expected electric field 
in an antenna (see, Table 1 from Murphy et al, 2008). 

• In case of the oblique magnetic field 9f )W = 45° the magnetic barrier is formed at 
the front of the heat shield. 

• The simulations were performed for the distance r = 4.5 R s because of a 
numerical limitation for convergence in the iterations of the electric field equation 
(Lipatov, 2002). For distance r = 9.5 R s one needs to design another algorithm for 
solving the electric field equation, Eq. 8. We are now trying to use the implicit 
’’splitting of operator” method to solve the generalized Ohm’s law equation, Eq. 8. 

• Our computations were performed only for the case (5 e = 0.00125 — 0.025 that is 
smaller at least by factor of 4 than realistic values. For higher value of f3 e one 
needs computation with significantly larger number of the macro-particles. It may 
be possible with the use of more than 320 processors on ” SGI-Columbia” parallel 
computer system (NASA Ames Center) with MPI environment to achieve a good 
scalability. 

• All simulations for 6b u = 11° were performed for a short period physical time 
because the present code needs a lot of time steps to compute the electric field. 

A more effective algorithm needs to be designed for longer period of the physical 
time. It will be important for an investigation of the low plasma density cavern 


formation around the SP+ bus. 
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• Simulated electric field perturbations are comparable to or exceed maximum 
electric field one expects for the SP+ spacecraft. Therefore, a wake design of SP+ 
electric field plasma-wave antenna may not be a viable option for SP+. 

Future work that is still needed: 
ffybrid simulations 

• Simulation with higher value of fi e — 0.1, higher value of the solar wind velocity 
U 0 ~ 400 — 800km/s, and lower value of the magnetic field B 0 « 1500 nT. 

• Simulations with various boundary condition on the spacecraft bus and trusses 
configuration; 

• Simulation with much higher space resolution and a number of macro-particles to 
reduces a ’’shot” noise in density and polarization electric field; 

• Longer simulations in order to receive a ’’steady-state” solution with a formation 
of a plasma wake in quasi-parallel interaction. 

The results of these simulations will be included in the final Report. 

Full kinetic simulations 

• Simulation that will take into account the charging of the spacecraft, the charge 
separation effects, an outgassing from heat shield, a photoionization and an 
electron impact ionization effects near the spacecraft. 
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Table 1 . Important values of the electromagnetic field in an antenna from 
” Thermal Design Consideration for the Solar Probe Plus Electric Field An- 
tenna by Murphy et al., Feb. 2008”. 



4 R s 

12i? s 

20 R s 

U xB 

12 V/m 

2.5 V/m 

0.5 V/m 

max i DC 

lV/m 

0.1 V/m 

0.1 V/m 

Emax, 100kHz 

10“ 3 V/m Hz” 1 / 4 

10" 4 V/m Hz" 1 / 4 

10" 4 V/m Hz" 1 / 4 

V Noise at 100 kHz 

10“ 6 V/Hz- 1 / 4 

3 x 10 -7 V/Hz -1 / 4 

2 x 10 -7 V /Hz -1 / 4 
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Table 2. Solar Wind Parameters for Modeling. 


Case Uq 

nsw 

^0 

@bu 

M a P p 

Pe 

a 

200 km/s 

10 4 

6250 nT 

n° 

1.2 

0.1 

0.000125 

b 

200 km/s 

10 4 

6250 nT 

n° 

1.2 

0.1 

0.025 

c 

200 km/s 

10 4 

6250 nT 

45° 

2.0 

0.1 

0.000125 
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Table 3. Electromagnetic Field Perturbation in Plasma Environment near SP + . 


Case 

t(T ce ) 

5E X (1.24: V/m) 

SE y ( 1.24 V/m) 

(5E 2 (1.24 V/m) 

SB X (6250 nT) 

SB y (6250 nT) 

SB z (6250 nT) 

E n at bus 

a 

0.014 

0.005-0.01 

0.01-0.5 

0.05-0.1 

0.001 

0.0002 

0.001 

1000 

a 

0.29 

0.005-0.01 

0.05-1.0 

0.05-0.2 

0.005 

0.01 

0.1 

1200 

b 

0.014 

0.002-0.01 

0.05 

0.05 

0.001 

0.0002 

0.001 

30 

b 

0.03 

0.02-0.04 

0.05 

0.02-0.2 

0.005 

0.01 

0.04 

30 

c 

13.6 

0.2-0. 5 

0.25-0.5 

0.2-2. 5 

0.05-1 

0.05-2 

0.05- 1.4 

5.0 











Figure 1 . Scheme of the interaction of the solar wind with SP+. The spiral magnetic 
field is inside the x-z plane. (1) a heat shield; (2) a support system; (3) a spacecraft bus. 
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Figure 2. Solar wind ion density in the y-x ( z = 0) and z-x (y = 0) planes. L 
200km/s, Ma = 1.5, B 0 = 6250 nT, E 0 = UqB 0 = 1.24V/m, nsw = 10 4 cm~ 3 , f3 v = 
P e = 0.000125, 9bu = 11°. Linear perturbations at the beginning of simulation 
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Figure 3. Electric field component E x ( 1.24 V/m) in the y-x (z = 0) and z-x (?/ = 0) 
planes. Uq = 200km/s, = 1.5, B 0 = 6250 nT, nsw — 10 4 cm -3 , /3 P = 0.1, /3 e = 

0.000125, 9b u = 11°. Linear perturbations at the beginning of simulation t = 0.014 T ce 
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Figure 4. 1-D cuts for the electric field component E x ( 1.24 V/m) (top: y = 0, z = 0; 
middle: x = 1.5, z = 0; bottom: x = 1.5, y = 0). U Q = 200km/s, M A = 1.5, 
Bo = 6250 nT, risw = 10 4 cm -3 , /3 P = 0.1, f3 e = 0.000125, 0} m = 11°. Linear perturbations 
at the beginning of simulation t = 0.014T ce 
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Figure 5. Electric field component E y ( 1.24 V/m) in the y-x (z = 0) and z-x (?/ = 0) 
planes. Uq = 200km/s, Ma = 1.5, i?o = 6250 nT, nsw = 10 4 crn~ 3 , /3 P = 0.1, /3 e = 
0.000125, 9b u = 11°. Linear perturbations at the beginning of simulation t = 0.014 T ce 
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Figure 6. 1-D cuts for the electric field component E y ( 1.24V/m) (top: y = 0, z = 0; 
middle: x = 1.5, z — 0; bottom: x = 1.5, y — 0). U 0 — 200km/s, M A = 1.5, 
Bo = 6250 nT, risw = 10 4 cm -3 , /3 P = 0.1, f3 e = 0.000125, 0} m = 11°. Linear perturbations 
at the beginning of simulation t = 0.014 T ce 
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Figure 7. Electric field component E Z {1.2A /m) in the y-x (z = 0) and z-x (y = 0) planes. 
U 0 = 200km/s, M A = 1.5, B 0 = 6250 nT, n sw = 10 4 cm~ 3 , /3 P = 0.1, (3 e = 0.000125, 
9b u = 11°. Linear perturbations at the beginning of simulation t = 0.014 T ce 
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Figure 8. 1-D cuts for the electric field component E~( 1.24 V/m) (top: y = 0, z = 0; 
middle: x = 1.5, z — 0; bottom: x = 1.5, y = 0). U 0 = 200km/s, M A = 1.5, 
Bo = 6250 nT, risw = 10 4 cm -3 , /3 P = 0.1, f3 e = 0.000125, 0} m = 11°. Linear perturbations 
at the beginning of simulation t = 0.014 T ce 
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Figure 9. Magnetic field component 5 x (6250nT) in the y-x ( z = 0) and z-x (y = 0) 
planes. Uq = 200km/s, Ma = 1.5, B 0 = 6250 nT, nsw = 10 4 cm -3 , /3 P = 0.1, /3 e = 
0.000125, 9b u = 11°. Linear perturbations at the beginning of simulation t = 0.014 T ce 
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Figure 10. Magnetic field component 1^(6250 nT) in the y-x {z = 0) and z-x ( y = 0) 
planes. Uq = 200km/s, Ma = 1.5, B 0 = 6250 nT, nsw — 10 4 cm -3 , /3 P = 0.1, /3 e = 
0.000125, 9b u = 11°. Linear perturbations at the beginning of simulation t = 0.014 T ce 
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Figure 11. Magnetic field component _B z (6250nT) in the y-x {z = 0) and z-x ( y = 0) 
planes. Uq = 200km/s, Ma = 1.5, B 0 = 6250 nT, nsw — 10 4 cm -3 , /3 P = 0.1, /3 e = 
0.000125, 9b u = 11°. Linear perturbations at the beginning of simulation t = 0.014 T ce 
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Figure 12. Solar wind ion density in the y-x (z = 0) and z-x ( y = 0) planes. 
Uq = 200km/s, Ma = 1.5, B 0 = 6250 nT, E 0 = UoB 0 = 1.24V/m, nsw = 10 4 cm -3 , 
/ 3 P = 0.1, /3 e = 0.000125, 6f m = 11°. Nonlinear saturation of the perturbations at 
t = 0.29T ce (0.04:T transit ) 
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Figure 13. Electric field component F7 x (1.24V/m) in the y-x (z = 0) and z-x 
(y = 0) planes. Uq = 200km/s, Ma = 1.5, B 0 = 6250 nT, nsw = 10 4 cm -3 , 
(3 P = 0.1, /3 e = 0.000125, 9f, u = 11°. Nonlinear saturation of the perturbations at 


t — 0.29 T ce ( 0 . 04 T transit ) 




Figure 14. 1-D cuts for the electric field component E x ( 1.24 V/m) (top: y = 0, z = 0; 
middle: x = 1.5, z — 0; bottom: x = 1.5, y — 0). U 0 — 200km/s, M A = 1.5, 
Bo = 6250 nT, n = 10 4 cm -3 , /3 P = 0.1, (3 e = 0.000125, 9b u = 11°. Nonlinear saturation 
of the perturbations at t = 0.29 T ce (0.04T transit ) 
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Figure 15. Electric field component F7 y (1.24 V/m) in the y-x {z = 0) and z-x 
(y = 0) planes. Uq = 200km/s, Ma = 1.5, B 0 = 6250 nT, nsw = 10 4 cm -3 , 
(3 P = 0.1, /3 e = 0.000125, 9f, u = 11°. Nonlinear saturation of the perturbations at 



t — 0.29 T ce ( 0 . 04 T transit ) 
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Figure 16. 1-D cuts for the electric field component E y ( 1.24V/m) (top: y = 0, z = 0; 
middle: x = 1.5, z — 0; bottom: x = 1.5, y — 0). U 0 — 200km/s, M A = 1.5, 
Bq = 6250 nT. risw = 10 4 cm -3 , (3 V = 0.1, /3 e = 0.000125, 9b u = 11°. Nonlinear saturation 
of the perturbations at t — 0.29 T ce (0.04T transit ) 
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Figure 17. Electric field component _E z (1.24V/m) in the y-x (z = 0) and z-x 
(y = 0) planes. Uq = 200km/s, Ma = 1.5, B 0 = 6250 nT, nsw = 10 4 cm -3 , 
(3 P = 0.1, /3 e = 0.000125, 9f, u = 11°. Nonlinear saturation of the perturbations at 


t — 0.29 T ce ( 0 . 04 T transit ) 
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Figure 18. 1-D cuts for the electric field component E~( 1.24 V/m) (top: y = 0, z = 0; 
middle: x = 1.5, z = 0; bottom: x = 1.5, y — 0). U 0 — 200km/s, M A = 1.5, 
Bq = 6250 nT, risw = 10 4 cm -3 , /3 P = 0.1, (5 e = 0.000125, Q\ )U = 11°. Nonlinear saturation 
of the perturbations at t — 0.29 T ce (0.04T transit ) 
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Figure 19. Magnetic field component 5 a ,(6250nT) in the y-x (z = 0) and z-x 
(y = 0) planes. Uq = 200km/s, Ma = 1.5, B 0 = 6250 nT, nsw = 10 4 cm -3 , 
(3 P = 0.1, {3 e = 0.000125, Obu = 11°. Nonlinear saturation of the perturbations at 


t — 0.29 T ce (0.04 T transit ) 
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Figure 20. Magnetic field component 5^(6250 nT) in the y-x (z = 0) and z-x 
(y = 0) planes. Uq = 200km/s, Ma = 1.5, B 0 = 6250 nT, nsw = 10 4 cm -3 , 
(3 P = 0.1, {3 e = 0.000125, 9f, u = 11°. Nonlinear saturation of the perturbations at 
t = 0.29 T ce ( 0 . 04 T transit ) 
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Figure 21. Magnetic field component 5,(6250 nT) in the y-x ( z = 0) and z-x 
(y = 0) planes. Uq = 200km/s, Ma = 1.5, B 0 = 6250 nT, nsw = 10 4 crif 3 , 
(3 P = 0.1, /3 e = 0.000125, 9f, u = 11°. Nonlinear saturation of the perturbations at 
t = 0.29 T ce ( 0 . 04 T transit ) 
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Figure 22. Solar wind ion density in the y-x (z = 0) and z-x ( y = 0) planes. U 0 = 
200km/s, Ma = 1-5, B 0 = 6250nT, E 0 = UqB 0 = 1.24V/m, nsw = 10 4 cm~ 3 , = 0.1, 

P e = 0.025, 9bu = 11°. Linear perturbations at the beginning of simulation t = 0.014 T ce 
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Figure 23. Electric field component E x ( 1.24 V/m) in the y-x (2 = 0) and z-x (y = 0) 
planes. W = 200km/s, Ma = 1.5, B 0 = 6250 nT, nsw = 10 4 cm -3 , (3 P = 0.1, (3 e = 0.025, 
dbu = 11°. Linear perturbations at the beginning of simulation t = 0.014 T ce 
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Figure 24. 1-D cuts for the electric field component E x { 1.24 V/m) (top: y = 0, z = 0; 
middle: a; = 1.5, z — 0; bottom: x = 1.5, y — 0). U 0 — 200km/s, M A = 1.5, 
Bq = 6250 nT. risw = 10 4 cm -3 , /3 P = 0.1, j3 e = 0.025, 6b u = 11°. Linear perturbations 
at the beginning of simulation t = 0.014 T ce 
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Figure 25. Electric field component £^(1.24 V/m) in the y-x (z = 0) and z-x (y = 0) 
planes. Uq = 200km/s, Ma = 1.5, i?o = 6250 nT, nsw = 10 4 cm -3 , (3 P = 0.1, f3 e = 0.025, 
9b u = 11°. Linear perturbations at the beginning of simulation t = 0.014 T ce 
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Figure 26. 1-D cuts for the electric field component E y ( 1.24 V/m) (top: y = 0, z = 0; 
middle: x = 1.5, z = 0; bottom: x = 1.5, y — 0). U 0 — 200km/s, M A = 1.5, 
B 0 = 6250 nT, = 10 4 cm- 3 , = 0.1, & = 0.025, = 11°. Linear perturbations 

at the beginning of simulation t = 0.014 T ce 
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Figure 27. Electric field component F z (1.24V/m) in the y-x (z = 0) and z-x (y = 0) 
planes. Uq = 200km/s, Ma = 1.5, B 0 = 6250 nT, nsw = 10 4 cnf 3 , (3 P = 0.1, (3 e = 0.025, 
dbu = 11°. Linear perturbations at the beginning of simulation t = 0.014 T ce 





Figure 28. 1-D cuts for the electric field component F 2 (1.24 V/m) (top: y = 0, z = 0; 
middle: x = 1.5, z — 0; bottom: x = 1.5, y — 0). U 0 — 200km/s, M A = 1.5, 
Bq = 6250 nT, risw = 10 4 cm -3 , /3 P = 0.1, j3 e = 0.025, 6b u = 11°. Linear perturbations 
at the beginning of simulation t = 0.014 T ce 
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Figure 29. Magnetic field component _B x (6250nT) in the y-x {z = 0) and z-x (y = 0) 
planes. Uq = 200km/s, Ma = 1.5, B 0 = 6250 nT, nsw = 10 4 cm~ 3 , (3 P = 0.1, f3 e = 0.025, 
dbu = 11°. Linear perturbations at the beginning of simulation t = 0.014 T ce 
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Figure 30. Magnetic field component 1^(6250 nT) in the y-x (z = 0) and z-x (y = 0) 
planes. Uq = 200km/s, Ma = 1.5, B 0 = 6250 nT, nsw = 10 4 cm -3 , /3 P = 0.1, (3 e = 0.025, 
dbu = 11°. Linear perturbations at the beginning of simulation t = 0.014 T ce 
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Figure 31. Magnetic field component _B z (6250nT) in the y-x {z = 0) and z-x (y = 0) 
planes. Uq = 200km/s, Ma = 1.5, B 0 = 6250 nT, nsw = 10 4 cm~ 3 , /3 P = 0.1, (3 e = 0.025, 
dbu = 11°. Linear perturbations at the beginning of simulation t = 0.014 T ce 
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Figure 32. Solar wind ion density in the y-x (z = 0) and z-x ( y = 0) planes. U 0 = 
200km/s, Ma = 1-5, B 0 = 6250nT, E 0 = UqB 0 = 1.24V/m, nsw = 10 4 cm~ 3 , = 0.1, 

P e = 0.025, 9bu = 11°. Nonlinear saturation of the perturbations at t — 0.03 T ce 
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Figure 33. Electric field component _E x (1.24V/m) in the y-x (z = 0) and z-x (y = 0) 
planes. Uq = 200km/s, Ma = 1.5, B 0 = 6250 nT, ngiy = 10 4 cm~ 3 , /3 P = 0.1, (3 e = 0.025, 
dbu = 11°. Nonlinear saturation of the perturbations at t — 0.03 T ce 
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Figure 34. 1-D cuts for the electric field component E x ( 1.24 V/m) (top: y = 0, z = 0; 
middle: x = 1.5, z — 0; bottom: x = 1.5, y — 0). U 0 — 200km/s, M A = 1.5, 
Bq = 6250 nT. risw = 10 4 cm -3 , j3 v = 0.1, /3 e = 0.025, 9b u = 11°. Nonlinear saturation of 
the perturbations at t — 0.03 T ce 
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Figure 35. Electric field component E y ( 1.24V/m) in the y-x (z = 0) and z-x (y = 0) 
planes. Uq = 200km/s, Ma = 1.5, B 0 = 6250 nT, nsw = 10 4 cm -3 , (3 P = 0.1, f3 e = 0.025, 
9b u = 11°. Nonlinear saturation of the perturbations at t — 0.03 T ce 
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Figure 36. 1-D cuts for the electric field component E y { 1.24V/m) (top: y = 0, z = 0; 
middle: x = 1.5, z — 0; bottom: x = 1.5, y — 0). U 0 — 200km/s, = 1.5, 

i?o = 6250 nT. risw = 10 4 cm -3 , (3 V = 0.1, /3 e = 0.025, 9b u = 11°. Nonlinear saturation of 
the perturbations at t — 0.03 T ce 
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Figure 37. Electric field component F z (1.24V/m) in the y-x (z = 0) and z-x ( y = 0) 
planes. Uq = 200km/s, Ma = 1.5, B 0 = 6250 nT, risw = 10 4 cm -3 , (3 P = 0.1, (3 e = 0.025, 
dbu = 11°. Nonlinear saturation of the perturbations at t — 0.03 T ce 
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Figure 38. 1-D cuts for the electric field component F 2 (1.24 V/m) (top: y = 0, z = 0; 


middle: x = 1.5, z — 0; bottom: x = 1.5, y — 0). U 0 — 200km/s, M A = 1.5, 


£>o = 6250 nT, risw = 10 4 cm 3 , (3 P = 0.1, /3 e = 0.025, 9b u = 11°. Nonlinear saturation of 


the perturbations at t = 0.03 T ce 
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Figure 39. Magnetic field component _B x (6250nT) in the y-x (z = 0) and z-x (y = 0) 
planes. Uq = 200km/s, Ma = 1.5, B 0 = 6250 nT, nsw = 10 4 cm -3 , (3 P = 0.1, f3 e = 0.025, 
dbu = 11°. Nonlinear saturation of the perturbations at t — 0.03 T ce 
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Figure 40. Magnetic field component 5^(6250 nT) in the y-x (z = 0) and z-x (y = 0) 
planes. Uq = 200km/s, Ma = 1.5, B 0 = 6250 nT, nsw = 10 4 cm -3 , (3 P = 0.1, (3 e = 0.025, 
dbu = 11°. Nonlinear saturation of the perturbations at t — 0.03 T ce 
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Figure 41. Magnetic field component £b(6250 nT) in the y-x (z = 0) and z-x (y = 0) 
planes. Uq = 200km/s, Ma = 1.5, B 0 = 6250 nT, riw = 10 4 cm -3 , /3 P = 0.1, (3 e = 0.025, 
dbu = 11°. Nonlinear saturation of the perturbations at t — 0.03 T ce 
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Figure 42. Solar wind ion density in the y-x (z = 0) and z-x ( y = 0) planes. U 0 = 
200km/s, M a = 2, B 0 = 6250 nT, E 0 = UqBq = 1.24V/m, nsw = 10 4 cm~ 3 , j3 v = 0.1, 
(3 e = 0.0125, 9 bu = 45°, t = 13.6 T ce (1.875 T transit . ) 
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Figure 43. Electric field component E x ( 1.24 V/m) in the y-x ( z = 0) and z-x (y = 0) 
planes. Uq = 200km/s, Ma = 2 , B 0 — 6250 nT, = 10 4 cm -3 , /? p = 0.1, (3 e = 0.0125, 
9 bu = 45°, t = 13.6 T ce (l. 875 T transit .) 
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Figure 44. Electric field component E y ( 1.24 V/m) in the y-x (z = 0) and z-x (?/ = 0) 
planes. Uq = 200km/s, Ma = 2 , B 0 = 6250 nT, nsw = 10 4 cm -3 , (3 P = 0.1, (3 e = 0.0125, 
9 bu = 45°, t = 13.6 T ce (l. 875 T transit .) 
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Figure 45. Electric field component E z ( 1.24 V/m) in the y-x ( z = 0) and z-x (y = 0) 
planes. Uq = 200km/s, Ma = 2 , B 0 — 6250 nT, = 10 4 cm -3 , /l p = 0.1, (3 e = 0.0125, 
9 bu = 45°, t = 13.6 T ce (l. 875 T transit .) 



Z/L Y/L 


Bx 




-2 0 2 4 

X/L 



-2 0 2 4 

X/L 

Figure 46. Magnetic field component B x in the y-x ( z = 0) and z-x (y = 0) planes. 
Uq = 200km/s, M a = 2, B 0 = 6250 nT, n sw = 10 4 cm- 3 , /3 P = 0.1, (5 e = 0.0125, 


Obu = 45°, t = 13.6 T ce (l. 875 T tranait .) 





Z/L Y/L 


4 


By 


0.04 



-2 0 2 4 

X/L 


By 



X/L 

Figure 47. Magnetic field component 5^(6250 nT) in the y-x (z = 0) and z-x (y = 0) 
planes. U$ = 200km/s, Ma = 2 , B 0 = 6250 nT, nsw = 10 4 cm -3 , f3 p = 0.1, /3 e = 0.0125, 
9 bu = 45°, t = 13.6 T ce (l. 875 T transit .) 
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Figure 48. Magnetic field component 5^(6250 nT) in the y-x (z = 0) and z-x ( y = 0) 
planes. U$ = 200km/s, Ma = 2, B 0 = 6250 nT, nsw = 10 4 cm -3 , f3 p = 0.1, /3 e = 0.0125, 
9 bu = 45°, t = 13.6 T ce (l. 875 T transit .) 
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